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LONG-TERM GOALS 

Our long-term technical goal is to produce a tested tangent linear and adjoint model for ROMS/TOMS 
(Regional Ocean Modeling System/Terrain-Following Ocean Modeling System) that is suitable for 
general use by ROMS/TOMS modelers. This is complementary to the Kalman Filter, ESSE, and 
Green’s Functions techniques being developed in other contexts. Our long-term scientific goal is to 
model and predict the mesoscale circulation and the ecosystem response to physical forcing in the 
various regions of the world ocean through ROMS/TOMS primitive equation modeling/assimilation. 

OBJECTIVES 

We seek to develop an adjoint model for ROMS/TOMS, a parallel/improved physics descendent of the 
serial SCRUM (Song and Haidvogel, 1994). We also seek to complete the assimilation system by 
including the adjoint in an estimation procedure for fitting the model to data. The resulting codes will 
be suitable for general use in any geometry of ROMS/TOMS, which presently lacks an adjoint. 

APPROACH 

This is fundamentally a collaborative effort involving University of Colorado (A. Moore), Rutgers (H. 
Arango) and Scripps (B. Comuelle, Ph. D. student E. Di Lorenzo, A. Miller, and D. Neilson). Our 
approach is to write the tangent linear and adjoint models for ROMS/TOMS by hand. With each 
participant in the project contributing expertise in coding and model testing, the approach is feasible. 
As the development is accomplished, the assimilation scheme is tested in various scenarios involving 
observations. The Scripps contingent will test the adjoint for ROMS/TOMS in the California Current 
CalCOFI region where they are presently applying ROMS (under NASA funding) to a physical- 
biological data synthesis and a model forecast scenario. Arango and Moore will test the adjoint in the 
Mid-Atlantic Bight (under NSF funding) for coupled atmosphere-ocean hindcast experiments using 
data collected at the observational network centered at the Long-Term Ecosystem Observatory (LEO- 
15). 

WORK COMPLETED 

The ROMS/TOMS adjoint team met many times over the past two years in intensive tangent linear and 
adjoint model writing/testing sessions. A working 2D and 3D tangent linear and adjoint model is now 
running and being used in various 2D and 3D applications. A large percentage of ROMS/TOMS 3D 
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routines are now completed and tested, but some of the less frequently used routines still need to be 
coded in for the adjoint. 

Finite-time normal modes have been tested in the 2D and 3D. The optimal perturbations scheme as 
been developed and tested in 2D and 3D. The stochastic optimals scheme has the driver completed but 
still needs testing. The 4D variational assimilation (4DVAR) is being constructed. So far, the driver 
for the inner and outer loops has been designed along with the modules needed to sample the model at 
observational locations and times. The cost function, pre-conditioning and descent algorithms are 
being built for 4DVAR. The ensemble prediction scheme needs the driver and module to compute 
random linear combinations of singular vectors. 

The 4DVAR and optimal perturbation schemes are now being set up to be tested in the Southern 
California Bight domain where CalCOFI hydrography and TOPEX altimetry observations are 
available. 

RESULTS 

The building of the adjoint is largely a technical task. Many scientific technical results will follow 
once the codes are ready. We here present an interesting application of the tangent linear and adjoint 
codes for a classic theoretical problem in physical oceanography. 

A 3D test case for computing optimal perturbations (Moore and Farrell, 1994) to study model error 
growth was set up in a simple stratified double gyre basin. The optimal perturbations are computed by 
an iterative application of forward and backward integrations of the tangent linear and adjoint 
equations in conjunction with the ARPACK library. Figure 1 shows an example from this case. The 
basic state surface circulation is shown in Figure la,b for a 10-day interval. This was generated by 
forcing the model for three years with steady sinusoidal winds and a relaxation to a hyperbolic tangent 
vertical temperature profile with linear latitudinal dependence. The circulation takes the form of two 
unsteady counter-rotating gyres that are often used as an analog of the subpolar/sub tropical gyre 
systems. The surface structure of the fastest growing optimal perturbation for this 10-day period is 
shown in Figure lc,d. Of all possible perturbations, this is the one that maximizes growth of 
perturbation energy over the 10-day period. At initial time, the sea level structure shows a pronounced 
upstream tilt against the basic state flow that is indicative of barotropic growth. As the perturbation 
evolves, it is advected downstream by the basic state circulation and grows in amplitude by a factor of 
40. 
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Figure 1: Panels a and b show the surface temperature and velocity vectors of the double gyre basic 
state flow 10 days apart. Panels c and d show the surface velocity and sea level at initial and final 
times for the optimal perturbation that maximizes the growth of perturbation energy over the 10 day 

period. 
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IMPACT/APPLICATIONS 


Users of the tangent linear and adjoint model for ROMS/TOMS will have a powerful tool for 
exploring data assimilation issues that include sensitivity to initial conditions and surface forcing, 
predictability and ocean dynamics. 

TRANSITIONS 

The work completed here will be part of the ROMS/TOMS utilities that will be freely available to all 
interested users. 

RELATED PROJECTS 

Moore, Arango, Miller and Cornuelle have a project funded by NSF (lead PI: A. Bennett, OSU) 
entitled “Modular Ocean Data Assimilation”. The goal is to use the infrastructure of the Inverse 
Ocean Modeling System of Chua and Bennett (2001) in conjunction with the ROMS/TOMS tangent 
linear and adjoint models for ocean data assimilation. The IOMS requires a somewhat different 
tangent linear model (but fortunately the same adjoint) that is based on the full fields rather than the 
perturbations fields. Both versions are now available for the 2D ROMS/TOMS kernel. Miller and 
Cornuelle are funded by NASA to explore Green’s Functions model fitting techniques (Miller and 
Cornuelle, 1999) with ROMS in the Southern California Bight of the California Current System 
http://osep.ucsd.edu/index.cgi7rsadjoint. Those results will prove useful in comparing with results 
from applying the adjoin t to these same data. 
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